On the speed of fluctuations around thermodynamic equilibrium 



Noah Linden, 1 Sandu Popescu, 2 Anthony J. Short, 3 and Andreas Winter 1,4 

1 Department of Mathematics, University of Bristol, University Walk, Bristol BS8 1TW, U.K. 
2 H.H.Wills Physics Laboratory, University of Bristol, Tyndall Avenue, Bristol BS8 1TL, U.K. 
Department of Applied Mathematics and Theoretical Physics, University of Cambridge, 
Centre for Mathematical Sciences, Wilberforce Road, Cambridge CB3 OWA, U.K. 
Centre for Quantum Technologies, National University of Singapore, 2 Science Drive 3, Singapore 117542 

(Dated:) 

We study the speed of fluctuation of a quantum system around its thermodynamic equilibrium 
state, and show that the speed will be extremely small for almost all times in typical thermodynamic 
cases. The setting considered here is that of a quantum system couples to a bath, both jointly 
described as a closed system. This setting, is the same as the one considered in [N. Linden et ai, 
Phys. Rev. E 79:061103 (2009)] and the "thermodynamic equilibrium state" refers to a situation 
that includes the usual thermodynamic equilibrium case, as well as far more general situations. 



Recently there has been significant progress in under- 
standing the fundamental principles of statistical me- 
chanics [J [2, H 0, H, H, 0, H ■ Underlying this progress is 
the realization that quantum mechanics allows individ- 
ual quantum states to exhibit statistical properties, and 
that ensemble or time averages are not needed to obtain 
a mixed equilibrium state for the system under consid- 
eration. This is a purely quantum phenomenon, and the 
key is entanglement, which leads to objective uncertainty 
- even when we have complete knowledge of the state of 
the whole system, a subsystem that is entangled with the 
rest of the system will be best described by a mixed state 
(i.e., a probabilistic mixture of pure states). 

This progress led to a proof from first principles that 
virtually any subsystem of any large enough system will 
reach equilibrium and fluctuate around it for almost all 
times [9j. In order to better understand the nature of 
these fluctuations, and to help understand the time-scales 
involved, here we investigate the speed of fluctuations 
around equilibrium: does the state of the subsystem os- 
cillate rapidly around equilibrium, or remain relatively 
static? Our main result is to put a universal upper bound 
on the average speed of fluctuations, showing that the 
speed is extremely small for almost all times in typical 
thermodynamic cases. 

Consider a large quantum system, described by a 
Hilbert space 7i. We decompose this system into two 
parts, a small subsystem S, and the rest of the system 
that we refer to as the bath B. Correspondingly, we 
decompose the total Hilbert space as TC = Hs ® Wb, 
where 7is and TLb (of dimensions ds and ds) are the 
Hilbert spaces of the subsystem and bath respectively. 
If the subsystem or bath would have infinite dimension, 
we bound its volume and energy to render the dimension 
finite, and project the interaction Hamiltonian onto the 
restricted Hilbert space H. 

Let the subsystem and bath evolve under a Hamilto- 
nian H that we decompose into a constant, subsystem, 



bath and interaction term 

H = H +H S + H B + H int (1) 

The decomposition is made unique by taking Hq propor- 
tional to the identity, Hs the tensor product of a traceless 
operator on the subsystem and the identity on the bath, 
Hb the tensor product of the identity on the subsystem 
and a traceless operator on the bath, and Hi nt trace- 
less for both the subsystem and bath. The total system, 
i.e. subsystem plus bath, is in a pure state |^(t)); let 
pit) = |\?(t))(\Ij(t)| be the density matrix representation 
of the state of the total system and let ps(t) = trg pit) 
be the density matrix of the subsystem. 

Following the notation in Q we define the time- 
averaged state of the whole system to, which is given by 

w = (p(t)) t = lim - [ p{t)dt. (2) 

Similarly we define ll>s = tr^w and u>b — trsa; as the 
time-averaged states of the subsystem and bath respec- 
tively. It is also convenient to introduce the notion of the 
effective dimension of a (mixed) state p by 



This is a measure of the number of states over which p is 
spread (e.g. for an equal mixture of N orthogonal states, 
d cfi = N). 

Clearly, wg is the canonical candidate for the equilib- 
rium state, but we need to pause to clarify just what it 
may mean that the system reaches equilibrium. Namely, 
note that due to the finiteness of the Hilbert spaces in- 
volved, there will be recurrences [on timescales exponen- 
tial incf ff (a;)], so a relaxation of ps(t) towards is out 
of the question. The best thing we can hope for in the 
current setting is that psit) remains close to los for most 
times t. 

To put our present results in a proper context, we first 
recall the key ideas and results of [9] . A key observation 
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was that the process of thermalisation actually contains 
many different aspects and we can decompose it into the 
following elements of analysis. 

1. Equilibration. We say that a system equilibrates if 
its state evolves towards some particular state (in gen- 
eral mixed) and remains in that state (or close to it) for 
almost all times. As far as equilibration is concerned, it 
is irrelevant what the equilibrium state actually is. 

2. Bath state independence. The equilibrium state of 
the system should not depend on the precise initial state 
of the bath, but only on its macroscopic parameters (e.g. 
its temperature) 

3. Subsystem state independence. If the subsystem is 
small compared to the bath, the equilibrium state of the 
subsystem should be independent of its initial state. 

4-. Boltzmann form of the equilibrium state. Under 
certain additional conditions on the Hamiltonian (espe- 
cially the interaction term) and on the initial state, the 
equilibrium state of the subsystem can be written in the 

familiar Boltzmannian form p$ — ^ exp ■ 

Realizing that thermalization can be decomposed in 
this way had major consequences. First, it allowed us 
to address each aspect separately. Second, and more im- 
portant, it allowed us to greatly expand the scope of our 
study. Indeed, we now consider equilibration as a gen- 
eral quantum phenomenon that may occur in situations 
other than those usually associated with thermalization. 
In particular we need not restrict ourselves to any of the 
following: standard thermal baths (that are described by 
a given temperature or restricted energy range), weak 
or short range interactions between the system and the 
bath, Boltzmannian distributions, situations in which en- 
ergy is an extensive quantity, etc. Furthermore, we can 
consider situations in which the subsystem does not reach 
equilibrium, and prove results about the bath or subsys- 
tem independence properties of the time-averaged state. 

In 0], we made substantial progress on items 1, 2, and 
3 above, under very weak assumptions. The only real 
constraint we impose on the Hamiltonian is that it has 
non-degenerate energy gaps. That is, any four energy 
eigenvalues satisfy E\ — E 2 = E$ — E4 if and only if 
Ei = E 2 and E 3 — E±, or E\ = E 3 and E 2 = E4. 
This assumption rules out non-interacting Hamiltonians 
of the form Hg 7^ 0, Hb ^ 0, H lnt — 0, which obviously 
do not equilibrate. These conditions are essentially those 
in Q but allow greater flexibility. Indeed, they allow the 
Hamiltonian to have degenerate energy levels, as long as 
the gaps between levels are non-degenerate, while in [9| 
degenerate levels were not allowed. Yet, it can be easily 
shown that even under these more general conditions the 
results of hold, namely that any small subsystem will 
reach equilibrium and fluctuate around it for almost all 
times [3-0( 1 . 

More precisely, Theorem 1] shows that the aver- 
age distance between ps(t) and its time average is 



bounded by 

^^•A{jk)A{^y < 4 > 

where D(p 1 ,p 2 ) = ~\\pi - p 2 \\i = \ tr yj (pi - p 2 ) 2 de- 
notes the trace distance between two density matrices 
(lo| . This is a natural distance measure on density ma- 
trices, giving the maximum difference in probability for 
any measurement on the two states . The meaning of 
eq. ([4]) is that the average distance between the instanta- 
neous state of the subsystem ps(t) and the fixed state ws 
will be small whenever the total dimension explored by 
the state (or the dimension explored in the bath) is much 
larger than the subsystem dimension. In typical thermo- 
dynamic situations, this will indeed be the case. Indeed, 
as dimensions typically grow exponentially with parti- 
cle number, we would expect any expression of the form 
poly(ds)/d c f{ (lu) to tend to zero in the thermodynamic 
limit (where the fraction of particles in the system tends 
to zero), as long as the energy distribution of |*(0)) was 
reasonable. This addresses point 1 in our programme; in 
practice one has to check that an initial state leads to 
large enough d cS (w); however, [ij Theorem 2] shows that 
this is the case for typical states from any sufficiently 
large subspace (for instance, a subspace of states with 
similar macroscopic properties) . 

Regarding item 2 it is proven 0, Theorem 3] that ini- 
tial states of tensor product form of a fixed state on the 
subsystem with a typical state of a large enough subspace 
on the bath, yield very similar equilibrium states. 

Finally, regarding item 3, we similarly prove 0, Theo- 
rem 3] that if the energy eigenstates of the Hamiltonian 
are sufficiently entangled, then also initial states of ten- 
sor product form of a typical state on the system with 
a fixed state on the bath have very similar equilibrium 
states. Furthermore, simple examples such as Hamiltoni- 
ans with tensor product eigenstates, show that without 
any additional assumptions on the correlation properties 
of the eigenstates one cannot expect subsystem indepen- 
dence to hold. 

We emphasize as in [§] that in the above discussion 
the "equilibrium" state ajg is not necessarily Boltzman- 
nian, and may depend on details of the Hamiltonian and 
the initial state (in particular, strong interactions, or con- 
served quantities for the subsystem will generate different 
equilibrium states). However, the equilibration results 
still hold. 

We now come to a crucial further issue, namely time 
scales. This issue needs to be a part of the general pro- 
gramme of investigating thermalisation. Again, this issue 
can be decomposed in a number of different questions. 
Firstly " How long does it take for a system to come close 
to equilibrium?", and secondly "What is the speed of 
fluctuations around equilibrium?" . This second point is 
what concerns us here. 
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While the magnitude of the fluctuations from equilib- 
rium may be small according to eq. ((4]) , this does not say 
anything about their speed. Mathematically, the speed 
of change of the subsystem state is given by 



vs(t) = lim 

(5t->0 



D(ps(t), Ps (t + 6t)) = 1 
St 2 



d Ps {t) 



dt 



(5) 



We will show that vs(t) is small for almost all times in 
typical thermodynamic cases, as follows. 

Theorem. For the average rate of change of ps it holds 
that 



(vs(t)) t < \\H s + H- mt \\< 



4 



d e s{u>) 



(6) 



where we take h= 1, and use the operator norm 11|. 



Proof. The time evolution of the subsystem state is given 
by 



dpsjt) 
dt 



= tr B (i[p(t),H])=J2 c k(t) 



(7) 



where e k with k = {1,2, d%} is a Hermitian orthonor- 
mal operator basis for the system, i.e. tr(efce^) = Ski- 
Hence, 

,.s , ( dpsjt) 
Ck (t) = tr s [-^e k 

= tt(i\p(t),H]e k ®I) 

= tr (p(t)i[H,e k ®I]) 

= tr(p(t)i[H s + H int ,e k ®I}). (8) 

Using our notation, and with a slight modification 
to use the operator norm, Reimann [7] shows that for 
a Hamiltonian with non-degenerate energy gaps, and a 
Hermitian observable A, 



Mp(t)A) - (tr(p(t)A)) t 



< 



\\Af 



t d eS {u)' 

Taking A = i{Hg + H- lntl e k ® I], and noting that 

\\i[Hs + H int ,e k ®I]\\ < 2\\H S + H int \\, 
we obtain 



, , . . n,2\ ^ M\Hs + I 

(c k {t) - (c k (t)} t ) ) < 



However, 



hence 



deff(w) 



(c fe (t)) t = (tr (t [p(t),H]e k (g>I)) t 
= tr (i [uj, H] e k ®I) = 0, 



{c k {t)Y) t < 



M\H S + H h 
d eS (uj) 



(9) 



(10) 



(11) 



(12) 



(13) 



This implies 



dp S (t) 



dt 



< 



^tr \^2c k (t)e k 

(e k ei) 

kl 

4|lg s + g int || 2 4 



(14) 



On the other hand, the trace norm and Hilbert-Schmidt 
norm are connected by the elementary relation \\X\\\ < 
(rankJC)||X|||. Hence, and using the concavity of the 
square root function, 



d Ps (t) 



dt 



a t 



< 



< 



\ 



dp s (t) 



dt 



\ 



ds 



d Ps {t) 



dt 



< 2\\H S + H U 



4 



d eS (uj)' 



(15) 



From the definition of vs(t) given by (O, we obtain the 
desired result 



(vs(t)) t = - 




(16) 



□ 



This result can be interpreted as follows. First, the 
speed of fluctuations varies in time and, of course, there 
may - and in general will - be times when the speed 
is extremely high. What our theorem shows is that, on 
average, the instantaneous speed is bounded by the ex- 
pression given in eq. ([6]) . Since speed is a positive quan- 
tity, this also means that the fraction of times for which 

v s (t) > K\\H S + Htatll^/ajfcy must be less than 1/K. 

Second, the speed of fluctuations depends linearly on 
the magnitude of the Hamiltonian, more precisely, on 
the subsystem and interaction Hamiltonian. Clearly 
the speed needs to depend linearly on the Hamiltonian 
since multiplying the Hamiltonian by a constant factor 
H — > AH makes the entire time evolution faster by the 
factor A. Furthermore, the speed of change of the state 
of the subsystem depends only on the part of the Hamil- 
tonian that acts directly on it, in particular it is inde- 
pendent of the bath Hamiltonian Hb- Of course, since 
the subsystem interacts with the bath, the time evolu- 
tion of the subsystem depends on the state of the bath 
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and thus implicitly on the bath Hamiltonian. However, 
as is already evident in eq. {SJ, the instantaneous change 
in the state of the subsystem (and hence the speed of its 
evolution) depends only on the instantaneous state of he 
bath and not directly on the bath Hamiltonian. Also, 
obviously, the time evolution is completely independent 
on the constant part of the Hamiltonian, Hq that only 
defines a reference point for the energy. These being said, 
the bound ([6|) should be better interpreted as a bound on 
the speed of fluctuations as measured in "natural units" 
i.e. 



H s + H 



< 



d% 



int 



(17) 



The main result of this paper is therefore that the av- 
erage speed, as measured in natural units, is bounded 

by yj d As mentioned earlier, because dimension 

generally grows exponentially with particle number, we 
would expect any fixed power of ds to be much smaller 
than d Q g(u>) in the thermodynamic limit. 

Note that similar arguments can be made for any fi- 
nite derivative of ps(t), exhibiting higher powers of the 
Hamiltonian in the upper bound with increasing degree 
of the time derivative. 

Finally, as far as the absolute value of the speed is 
concerned, we note that in most natural systems the 
magnitude of the Hamiltonian governing the speed, i.e. 
||.ffs + .Hint 1 1 scales relatively slowly (i.e. polynomially) 
with the number of particles. For example, in a system of 
n particles in which the Hamiltonian only contains two- 
particle interactions, we would expect the norm of the 
Hamiltonian to grow at most quadratically in n. Hence 
in the thermodynamic limit, when the total number of 
particles in the system increases, we would expect the ex- 
ponential dependence of the dimensional term \J d 
to dominate, causing the absolute value of the average 
speed to tend to zero. 

To summarize, we have shown that in the thermody- 
namic limit not only every subsystem spends almost all 
its time fluctuating very closely around a fixed state - 
the equilibrium state - but also that the speed of fluctu- 
ations becomes vanishingly small. Both these results at 
first glance appear to contradict the established view that 
a system has non-vanishing fluctuations around equilib- 
rium. However, as we have already argued in [9(, the 
global state ps(t) of a subsystem fluctuates extremely 
little around the equilibrium state, and what are gener- 
ally thought to be time fluctuations are just the result of 



quantum uncertainty (the probabilistic nature of the out- 
comes of quantum measurements) , that would be present 
even if the state were absolutely time independent. 

As in [9], our results hold not only for the standard 
statistical mechanical setting with Boltzmannian equi- 
librium, but under extremely general conditions. 
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